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PAST  A BACKWARD  FACING  STEP 

Isothermal  flow  and  flows  with  heat  transfer  and  property  variations 


RAVIKANTH  AVANCHA  AND  RICHARD  PLETCHER 
Department  of  Mechanical  Engineering 
Iowa  State  University,  Ames  I A 50011 


Abstract 

The  heat  transfer  and  fluid  mechanics  of  a turbulent  separating  and  reat- 
taching flow  past  a single-sided  backward-facing  step  are  studied  using 
large  eddy  simulation.  A fully  coupled,  low  Mach  number  preconditioned, 
collocated-grid,  central  differenced,  compressible,  finite  volume  formulation 
was  developed  to  conduct  the  simulations.  A sixth-order  compact  filter  was 
used  to  prevent  pressure- velocity  decoupling.  A compressible  version  of  the 
dynamic  subgrid  scale  model  was  used  to  model  the  effects  of  the  smaller  ed- 
dies. Navier-Stokes  characteristic  boundary  conditions  designed  by  Poinsot 
and  Lele  were  used  to  provide  boundary  conditions.  The  isothermal  tur- 
bulent flow  past  the  step,  at  a Reynolds  number  of  5,540  (based  on  the 
step  height  and  upstream  centerline  velocity)  and  a Mach  number  of  0.006, 
was  simulated  to  validate  the  formulation.  Subsequently,  the  bottom  wall 
downstream  of  the  step  was  supplied  with  constant  wall  heat  flux  levels 
of  1.0,  2.0,  and  3.0  kW/m2.  The  viscous  sub-layer  played  a critical  role 
in  controlling  the  heat  transfer  rate.  Streamwise  and  wall-normal  turbulent 
heat  fluxes  were  of  the  same  order  of  magnitude.  The  Reynolds  analogy  did 
not  hold  in  the  recirculation  region.  However,  the  Stanton  number  profiles 
showed  a striking  similarity  with  the  fluctuating  skin-friction  profiles. 

1.  Introduction 

Turbulent  flows  with  separation  and  reattachment,  in  the  presence  of  heat 
transfer,  occur  routinely  in  aircraft  propulsion  equipment,  such  as  turbines 
and  combustors,  and  cause  large  variations  of  local  heat  transfer  coefficient 
as  well  as  augmentation  of  overall  heat  transfer.  The  backward  facing  step 
geometry  is  well  suited  for  investigation  into  the  characteristics  of  flow  sep- 
aration, reattachment  and  recovery  in  the  presence  of  an  adverse  pressure 
gradient.  Most  of  the  numerical  calculations  to  date  for  this  flow  regime  with 
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heat  transfer  have  used  two  equation  turbulence  models  like  the  k — e and 
k — u models  with  the  Reynolds  averaged  Navier-Stokes  (RANS)  equations. 
Existing  RANS  approaches  have  not  been  successful  in  predicting  all  the 
flow  features  and  heat  transport  mechanisms  for  this  geometry.  A summary 
of  benchmark  calculations  [1]  noted  that  only  the  results  for  streamwise 
velocity  and  dissipation  agreed  amongst  all  contributors,  but  substantial 
differences  were  noted  for  the  wall-normal  velocity,  skin  friction,  Nusselt 
number,  and  wall  temperatures.  For  the  Nusselt  number  profiles,  none  of 
the  methods  produced  a shape  similar  to  the  experimental  results  of  [2]. 
There  did  not  appear  to  be  any  consistency  in  the  computed  points  of  max- 
imum Nusselt  number  and  the  reattachment  point  [1].  Many  simulations  in 
this  study  did  not  include  the  effects  of  property  variations  although  large 
heat  flux  levels  were  used.  Since  RANS  methods  rely  to  a great  extent 
on  the  modeling  aspects,  large  variations  in  the  quality  of  simulations  are 
observed  depending  on  the  model  assumptions.  Little  ad-hoc  modeling  is 
employed  in  large  eddy  simulation  (LES),  where  the  three-dimensional  un- 
steady motion  of  the  larger  eddies  is  computed  and  subgrid-scale  modeling 
is  employed  to  account  for  the  effects  of  the  smaller  eddies.  Thus,  the  large 
eddy  simulations  presented  in  this  work  offer  a way  to  better  understand 
the  effects  of  the  heat  transfer  and  property  variations  in  separated  flow 
regions  at  low  speeds. 


2.  Governing  Equations  and  Subgrid  Scale  Modeling 

The  Favre-filtered  compressible  Navier-Stokes  (NS)  equations,  and  the  equa- 


tion  of  state,  in  their  non-dimensional  form  are  given  as: 

dp  t dpUj 
dt  dxj 

= 0 

(1) 

dpUi  dpUiUj 

dp 

doij 

drij 

(2) 

dt  * dxj 

dxi 

dxj 

dxj 

+ ±o*& 

II 

m 

i 

dQj 

dxj 

(3) 

P 

= pRT 

(4) 

where  the  viscous  stress  tensor  and  heat  conduction  vector  are  given  as, 
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The  turbulent  stress  Txj  and  the  turbulent  heat  flux  Qj  have  to  be  modeled 
in  order  to  close  the  system  of  equations.  The  compressible  extension  [3]  to 
the  dynamic  subgrid  scale  model  [4]  employing  the  Smagorinsky  model  as 
the  base  model  has  been  used  in  this  study.  The  compressible  flow  version 
of  the  Smagorinsky  model  is  given  as 

Tij  = ^Tkkfiij  — 2 iiT{Sij  — -Skkdij)  ; j>%  = 2 ( + ~dxi)  } ^ 

where  pT  is  the  eddy-viscosity,  and  Sij  is  the  Favre-filtered  strain  rate  ten- 
sor. Tfcfc  = 2CipA2\S\2  is  used  to  model  the  sub-grid  scale  turbulent  kinetic 
energy,  r^,  as  proposed  by  [5].  For  closure,  pT  is  parameterized  by  equat- 
ing the  sub-grid  scale  energy  production  and  dissipation,  and  the  subgrid 
scale  heat  flux  vector  Qj  is  modeled  using  a gradient-diffusion  hypothesis, 
to  obtain 


„ pcsA2\s\df 

pT  = CspA  \J2bijbij  ; Qj  = — (8) 

where  Cs  is  a model  parameter  and  PrT  is  the  turbulent  Prandtl  number 
defined  as  the  ratio  of  eddy-viscosity  uT  to  eddy-diffusivity  aT.  Both  Cs  and 
PrT  are  computed  “dynamically”.  A is  the  filter  width  which  is  typically 
assumed  to  be  a function  of  the  grid  resolution,  and  calculated  as  Aav  = 
(AxAyAzY^ . 

3.  Numerical  Procedure 

A coupled  finite  volume  procedure,  in  primitive  variables  [p , u , v , w , T] 
was  used  to  solve  the  filtered  NS  equations.  The  method  was  fully  implicit, 
second  order  accurate  in  time,  with  advective  terms  disretized  using  second 
order  central  differences  and  viscous  terms  with  fourth  order  central  dif- 
ferences. Time  derivative  preconditioning  [6]  was  incorporated  to  alleviate 
the  stiffness  and  convergence  problems  associated  with  the  computation 
of  low  Mach  number  flows  using  traditional  compressible  formulations.  An 
all-speed  strategy  has  thus  evolved  that  enables  the  application  of  the  same 
methodology  to  incompressible  flows  and  compressible  flows  at  low  Mach 
numbers  where  effects  of  property  variations  need  to  be  accounted  for.  Sixth 
order  compact  filtering  [7]  was  used  to  eliminate  the  pressure-velocity  de- 
coupling peculiar  to  collocated-grid  methods.  The  system  of  algebraic  equa- 
tions was  solved  using  the  strongly  implicit  procedure  (SIP)  [8,  9].  A code 
well  optimized  for  performance  on  the  CRAY  T-90  was  used  to  perform  the 
large  eddy  simulations. 
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4.  Simulation  Details  and  Results 

The  computational  domain  for  the  large  eddy  simulations  was  designed  to 
match  the  backward-facing  step  geometry  from  the  particle  tracking  ve- 
locimetry  (PTV)  experiments  [10]  and  is  shown  in  Fig.l.  The  Reynolds 
number,  based  on  the  step  height  and  upstream  centerline  velocity  (which 
was  also  the  reference  velocity)  was  5,540.  Reference  values  of  thermal  con- 
ductivity, density  and  dynamic  viscosity  were  obtained  at  the  reference 
temperature,  Tref,  of  293  K.  The  grid  resolution  used  for  the  simulations 
was:  (1)  Upstream  of  step:  17  x 31  x 48,  (2)  Downstream  of  step:  72  x 46  x 48 
in  the  streamwise,  wall-normal  and  spanwise  directions  respectively.  Non- 
uniform  grids  were  employed  in  the  wall- normal  (y),  and  streamwise  ( x ) di- 
rections as  opposed  to  a uniform  grid  in  the  spanwise  (z)  direction.  No-slip 
boundary  conditions  were  enforced  at  the  top  and  bottom  solid  walls.  Pe- 
riodicity of  flow  was  assumed  in  the  spanwise  direction.  The  Navier-Stokes 
characteristic  boundary  condition  strategy  [11]  was  employed  at  the  inflow 
and  outflow  boundaries.  Turbulent  inflow  conditions  for  each  time  step  of 
the  simulation  were  provided  by  planes  of  data  stored  from  an  independent 
LES  of  a channel  flow  with  the  same  Reynolds  number  and  time  step.  For 
the  simulations  with  heat  transfer,  the  bottom  wall  downstream  of  the  step 
was  the  only  one  supplied  with  a constant  heat  flux.  The  remaining  walls 
were  insulated  (adiabatic  conditions). 

4.1.  ISOTHERMAL  FLOW  CASE 

The  streamwise  (Fig.  3),  wall- normal  (Fig.  4)  and  spanwise  mean  velocity 
distributions,  and  the  respective  root  mean  square  fluctuations  (Figs.  5,  6) 
from  the  simulation  showed  excellent  agreement  with  experimental  results. 
Third  order  moments  also  showed  good  qualitative  agreement  with  the  ex- 
periment. The  mean  reattachment  length  from  the  simulation  was  predicted 
to  be  6.0  x/h,  as  compared  with  6.51  x/h  reported  in  the  experiment.  The 
simulations  captured  the  presence  of  a smaller  corner  eddy  in  addition  to 
the  large  recirculation  bubble  (Fig.  2).  Details  of  the  isothermal  simulation 
and  results  have  been  presented  in  [12]. 

4.2.  CONSTANT  WALL  HEAT  FLUX  CASES 

Heat  flux  ( qw ) levels  of  1.0,  2.0,  and  3.0  kW/m2  corresponding  to  normal- 
ized heat  flux,  Q+  (=  qw/ pTefUTefCpTref)  levels  of  0.0014,  0.0028,  0.0042 
were  supplied  to  the  surface  downstream  of  the  step,  yielding  maximum 
Twaii/Tbuik  ratios  of  about  1.7,  2.3,  and  2.9,  respectively.  The  bulk  tem- 
perature profiles  (Fig.  7)  were  in  good  agreement  with  analytical  estimates 
based  on  an  overall  energy  balance  for  the  uniform  heat  flux  condition.  The 
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wall  temperatures  (Fig.  8)  showed  a dramatic  increase  downstream  of  the 
step,  and  reached  their  peak  values  in  the  neighborhood  of  the  streamwise 
distance  of  x/h  ~ 2.  This  increase  in  temperature  is  accompanied  by  a 
decrease  in  convective  heat  transfer  (as  evidenced  by  the  Stanton  number 
profiles  in  Fig.  9)  and  suggests  that  the  air  in  this  zone  is  almost  in  a “stag- 
nant” state.  The  Stanton  number,  St  (=  hj prefUrefCp)  is  a modified  Nus- 
selt  number  (St  = Nu/RePr ),  where  Nu  = hLy/kb,  Ly  is  the  step  height, 
h is  the  convective  heat  transfer  coefficient,  and  kb  is  the  bulk  thermal 
conductivity.  The  Stanton  number  attains  a maximum  slightly  upstream  of 
reattachment,  which  is  in  agreement  with  [13],  as  opposed  to  results  from 
several  other  studies  that  have  indicated  the  location  of  the  peak  Stan- 
ton (or  Nusselt  number)  to  coincide  with  the  mean  reattachment  point.  In 
the  region  of  reattachment,  the  impinging  shear  layer  is  responsible  for  the 
depression  in  the  wall  temperature  around  reattachment.  Downstream  of 
reattachment,  linear  increase  of  the  wall  temperature  and  the  monotonic 
decrease  in  the  mean  Stanton  number  profiles  is  consistent  with  the  growth 
of  the  thermal  boundary  layer  following  reattachment. 

The  Reynolds  analogy  does  not  hold  in  the  mean  sense,  i.e.,  the  mean 
Stanton  number  profiles  (Fig.  9)  do  not  correlate  with  the  mean  skin  friction 
profiles  (Fig.  10).  However,  it  is  interesting  to  note  that  the  mean  Stanton 
number  profiles  show  a more  striking  similarity  with  the  fluctuating  skin- 
friction  profiles  (Fig.  11)  than  they  do  with  the  profiles  of  the  average 
absolute  skin- friction  (Fig.  12),  or  the  mean  skin  friction.  The  high  degree 
of  correlation  between  the  fluctuating  skin-friction  coefficient  and  Stanton 
number  suggests  that  the  correct  velocity  scale  governing  the  strength  of 
convective  effects  in  the  reattachinent  zone  must  be  related  to  the  velocity 
fluctuations  rather  than  the  mean  velocity  [13]. 

The  bottom  wall  skin-friction  coefficient  (Fig.  10),  in  terms  of  normal- 
ized quantities  can  be  written  as  Cf  = (2 / Re)  [(n / p)  du / dy\w.  The  influ- 
ence of  the  wall  temperatures  for  the  three  heat  flux  cases  on  the  skin- 
friction  coefficient  is  through  the  density  and  viscosity  evaluated  at  the 
wall. 

It  is  indicated  in  [13]  that  the  streamwise  turbulent  heat  flux  is  negligible 
as  compared  to  the  wall-normal  turbulent  heat  flux.  However,  we  show  that 
the  streamwise  turbulent  heat  flux  (Fig.  13),  and  the  wall-normal  turbulent 
heat  flux  (Fig.  14)  are  of  the  same  order  of  magnitude.  Probable  reasons  for 
their  observation  of  negligible  streamwise  heat  flux  are  that  the  heat  flux 
levels  in  their  work  are  an  order  of  magnitude  smaller  and  the  Reynolds 
numbers  are  three  to  five  times  greater  than  the  ones  in  this  study. 

Root  mean  square  density  fluctuations  of  up  to  20  % were  observed.  At 
the  peak  wall  temperatures  for  the  three  heat  flux  cases,  the  viscosity  and 
thermal  conductivity  are  roughly  about  20%,  40%  and  60%  greater  than 
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the  values  at  the  inlet,  which  are  at  the  reference  temperature  of  293  K, 
thus  demonstrating  the  need  for  property  variations  to  be  considered  in 
calculations  involving  high  heat  fluxes. 

5.  Concluding  Remarks 

Large  eddy  simulations  to  study  the  heat  transfer  and  fluid  dynamics  of 
the  turbulent  reattaching  flow  past  a backward-facing  step  have  been  suc- 
cessfully conducted.  The  choice  of  the  formulation  enabled  the  inclusion 
of  property  variations,  and  facilitated  the  study  at  low  Mach  numbers. 
An  increase  in  heat  flux  supplied  to  the  bottom  wall  downstream  of  the 
step  results  in  the  heat  transfer  rate  starting  to  be  dominated  by  conduc- 
tion as  opposed  to  convection.  Large  variation  in  the  wall  temperature  is 
observed,  and  the  wall  temperature  attains  a maximum  close  to  the  step 
face.  The  peak  heat  transfer  rate  occurs  slightly  upstream  of  reattachment. 
Streamwise  and  wall-normal  turbulent  heat  fluxes  are  of  the  same  order  of 
magnitude.  The  Stanton  number  profiles  correlate  strongly  with  the  fluctu- 
ating skin-friction  profiles,  and  this  observation  underscores  the  importance 
of  the  near-wall  region  in  determining  the  heat  transfer  rate. 
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Figure  IS.  Streamwise  turbulent  heat  flux  Figure  14-  Vertical  turbulent  heat  flux 


